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AN ANALYSIS OF STATICALLY INDETERMINATE TRUSSES HAVING 
MEMBERS STRESSED BEYOND THE PROPORTIONAL LIMIT 
By Thomas W. Wilder, III 


SUMMARY 


A procedure for analyzing statically indeterminate trusses in the 
plastic stress range is presented vhich is applicable to trusses having 
any number of redundant members. A numerical example is used to i 'i 1 n r- 
trate the procedure. : - 


INTRODUCTION 


Although the truss has generally been replaced by the stiffene d 
sheet as the preferred structural form for the wing and fuselage, it — f — . “ 
still constitutes a useful configuration that finds application in air- 
craft-. The elastic stress analysis of the truss with redundant members 
is a classical problem with which all structural engineers are familiar. 
Aeronautical engineers, however, are interested in the maximum strength 
of trusses and therefore a plastic analysis of the structure is necessary. 
Two somewhat different procedures for analyzing statically indeterminate 
trusses with members stressed beyond the proportional limit have recently 
appeared in the aeronautical literature (refs. 1 and 2). Reference 1 
describes a procedure which In essence applies the principle of minimum 
complementary energy to plastic truss analysis but the method is limited 
to a single redundancy. The method of reference- 2 is an extension of the 
relaxation technique of R. V. Southwell and, In general, requires the 
determination of more unknowns than the static redundancy of the truss. 

Presented herein is an application of the principle of minimum com- 
plementary energy to the analysis of statically indeterminate trusses 
having members stressed beyond the proportional limit and having any num- 
ber of redundant members. The number of unknown quantities to be deter- 
mined is equal to the number of static redundancies. The use of the 
Ramberg and Osgood (ref. 3 ) analytical representation of the stress- strain 
curve that applies to aluminum alloys , magnesium, and stainless steel per- 
mits a concise analytic formulation of the problem for 'trusses made "of 
these materials. ... . 
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THEORY 


Of the possible combinations of forces, in the members of a truss, 
which satisfy equilibrium among themselves and are in equilibrium with 
the external loads, the correct ' combination is the one that minimizes the 
complementary energy (see ref. 4 for a discussion of the method of com- 
plement ary energy) : 

JE_ p^k , , 

u = e k (1) 

k=l ^0 


where 

U complementary energy 

m total number of members' in ‘truss 
e change' in length of a member 


P force acting in the member 

In addition to the usual assumptions of truss analysis, such as the 
assumption that the deformations are not sufficient to disturb the geo- 
metric relations between the members, the following assumptions are made: 

(1) Any movement of a support is perpendicular to the reaction at 
the support. 


(2) No forces are acting in the members before the application of 
the external loads. 

(3) The final stresses in the members are on the stress-strain 
curve and therefore have not resulted from plastic stress reversal. 


The integration indicated in equation (l) can be performed analytically 
through use of-the Ramberg-Osgood representation of the stress-strain 
curve (ref. 3)- This representation yields the following equation 
relating e and P: 


e 


= PL + 2 . P R L /F \ n 
AE 7 AE \P R/ ) 


( 2 ) 


where 

E elastic modulus 

L length of the member 


A 


cross-sectional area of the member 
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i P R force causing a stress in the member equal to the 0.7E secant yield 

stress 

* n shape parameter used by Kamberg and Osgood for the stress-strain 

curve . 

Relation (2) expresses the change in length as the sum of two terms: 
the first is recognized as the change in length that would occur if the 
member were elastic; the second term is the additional change in length 
caused by plastic deformation. 

Substituting equation (2) into equation (l) and subsequently inte- 
grating yields 


Let 


¥ 


U = 



pA + 3 fi c \ nk ~ 

^k^k 7 ( n k + 1 ) A k E k\PE k j 


(3) 


p k = p 0k + a l p lk + • • • + a*?^ + . - . + (4) 


where 

P Qk statically determinate force in the kth member caused 

by the external loading when J redundant members are 
assumed missing - 

P 1k statically determinate force in the kth member caused 

by a unit tensile force in the ith redundant member 
when the applied loads and the forces in the other 
redundant members are assumed equal to zero 


8.^^ • • ■ 8 j 


coefficients to be determined 


When equation (4) is substituted into equation (3) and the resulting 
expression for the complementary energy is minimized by equating to zero 
itB partial derivatives with respect to each a, there results a s et 
of j nonlinear simultaneous equations which determine a^, . . . aj. 
The general equation obtained by setting 



(i = 1, 2 , 


j) 


(5) 



k — 
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is 


^Lk^ik-^k , v m ~ p ik^k Pjk E ik E k 

/ + . . .+&js + • . . + a j / — * 

fa A k E k fa A k®k J fa A k E k 


jr- Pok^ikLk 

k=l Ajs^k 


3 y^~ ^^iA/Ppk + a lPjk + • ‘ • + ajPjk + • • • + ajPjiA ^ 

7 k=l A k E k \ P R k / 


( 6 ) 


(i * 2^^ . . . j) 


With the exception of the second summation on the right-hand side, equa- 
tion (6) is the same as would have resulted from an analysis of an elastic 
truss. The additional summation in equation (6) accounts for the effects 
due to plasticity. 


SOLUTION- OF THE NONLINEAR SIMULTANEOUS EQUATIONS 


A description of- a convenient procedure for solving the nonlinear 
simultaneous equations follows. (Experienced computers may prefer a 
method with which they are more familiar. ) 

The set of equations (6) may he written 


M'M = | k e| + l F ( a >i 


(7) 


|a| = [a]~ 1 |k e | + [A]" 1 |F(a)| 



or 
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where 

M 



square matrix having as elements the summations appearing as 
coefficients of the a's in the left-hand, sides of the 
equations 

inverse of [a] 

column matrix of the a’s 

column matrix consisting of the first (elastic) te rms in the 
right-hand sides of the equations 

column matrix of the second (plastic) terms in the right- hand 
sides of the equations 


Were stress and strain related by the elastic modulus for all forces 
in all members, the solution would he " 


a 




( 9 ) 


Equation (9) represents the elastic solution and may he used as a 
first approxi m ation to the values for the a'B which are to satisfy 
the nonlinear equations. The second approximation will then he 


|»n>| - |a<°>| + [a]-i|f(J°))| 

and the rth approximation after the first will he 

|a< r >| - |a<°)| ♦ [ A ]- 1 |F(a( r - 1 >)| 


( 10 ) 


( 11 ) 


Each set of values for the a 1 s should he compared with the pre- 
ceding values to determine whether or not continuance of the iterative 
procedure indicated by equation ( 11 ) is Justified. For trusses having 
several redundant members, repeated application of equation (11) may 
yield varying rates of convergence or divergence for successive values, 
of an a. If the iteration appears unsatisfactory, the procedure can 
usually he altered to produce rapid convergence. _ 
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NUMERICAL EXAMPLE 


For illustration of the procedure, a truss having two redundant 
members will be considered. If the stress-strain curve for the material 
of each member is the same, and if the tensile and compressive portions 
of the curve are symmetrical about the origin, equation (6) yields: 


m 


k=l 


F ih 2;L k + a p ik F 2k L k = _ y— WlA _ 

A 


k=l 


Av 


k=l 


F R k P lk L k /P Qk + + agPgkV 


T k=l 


> ( 12 ) 


m 


k=l 


F lk p 2k L k TT” p 2k 2l 'k _ 

+ ao / = “ 

A k k=l ' A k 


m 


F 0k p 2k L k 


k=l A k 


3 FR k F2kLk / p Ok + a l F lk + & 2 P 2k\ n 


T k=l A k 


■Ht 


For the statically indeterminate trass given In figure 1(a), the 
statically determinate load system is shown in figure 1(b) and the two 
redundant load systems are shown in figures 1(c) and 1(d). Lengths and 
areas of the truss members are summarized in table— I, along with the 
statically determinate forces. 

The stress-strain curve used was arbitrarily chosen and has a 0.7E 
secant yield stress of- 40.5 kips per square inch and a value of 6.56 
for n. 

The sums obtained from table II are used (as indicated below the 
table) to determine the elastic solution, which’ Is the first approxima- 
tion. The sums required for the second approximation are given in 
table III, and the application of . equation (10) to obtain the second 
approximation is indicated below table III. 
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Further iterations indicated an oscillatory divergence of the coef- 
ficients a-j_ and an oscillatory convergence of &£. Convergence was 

obtained when averages of successive values of the coefficient were 
used in the iteration. The final values for a^ and ag are used in 

table IV; convergence is indicated when successive solutions of eq.ua— 
tion (11) result in identical values. - - ' 


CONCLUDING REMABKS 


The analytical expression for the stress-strain curve used in the 
present method for analyzing statically indeterminate trusses having 
members stressed above the proportional limit has a form which permits' 
the analysis to be divided into two parts. The first part, or elastic 
solution, is simply an application of the familiar principle of least 
work to an elastic truss. The second part of the procedure is a routine 
iterative process which accounts for the plastic deformations of the 
truss and yields a correction to the elastic solution. 

The method is applicable to a truss having any number of redundant 
members. The number of unknown quantities to be determined is equal to - ” 
the number of redundant members in the truss, and the evaluation of these 
unknown quantities determines directly the force acting in each member 
of the truss. When the forces are known, the deformations can readily 
be found. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , September 16, 1952. 
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TABLE I.- DIMENSIONS AND STATICALLY DETERMINATE 
FORCES USED IN NUMERICAL EXAMPLE 


Member 

Length 

Area 

p o 

p l- 

p 2 

a 

30 ' 

0.25 

7-5 

- 0.6 

0 

b 

50 

.20 

12.5 

1.0 

0 

c 

50 

.20 

0 

1.0 

0 

d 

30 

.25 

- 15.0 

- 0.6 

0 

e 

4o 

.20 

0 

- 0.8 

- 0.8 

f 

30 

.25 

7.5 

0 

-0.6 

g. 

50 

.20 

0 

0 

1.0 . 

h 

50 

.20 

- 12.5 

0 

1.0 

I 

30 

.25 

0 

0 

>- 0.6 

2 

4o 

.20 

10.0 

0 

-0.8 
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TABLE II.- COMPUTATIONS FOR ELASTIC SOLUTION 
(FIRST APPROXIMATION) 


Member 

Pl2L 

PlP 2 L 

P 2 2 L : 



A 

A 

A 

warn 

A 

a 

43.2 

0 

0 

-540 

0 

b 

250.0 

0 

0 

3125 

0 

c 

250.0 

0 

0 

0 

0 

d 

43.2 

0 

0 

1080 

0 

e 

128.0 

128.0 

128.0 

0 

0 

f 

0 

0 

43.2 

0 

-54o 

g 

0 

0 

250.0 

0 

0 

h 

0 

0 

250.0 

0 

-3125 

i 

0 

0 

43.2 , 

0 

0 

J 

0 

0 

128.0 

0 

-1600 

z 

714.-4 

128.0 

842.4 

- 

3665 

-5265 


M - 1 

= ■ 

714.4 

128.0 

128.0 

824.4 

-1 

0.00143895 

0.00021864 

-0.00021864 

0.00122031 

|a<°> 

- [A] - 1 

-3665 





5265 


a]L (0> 


-6.4249 

a 2 (°) 


7. 2262 

















a 1 P> 


TABLE III.- C (MUTATIONS FOR SECOND APPROXIMATION 


Member 



10.125 
8.100 
- 8.100 
-10.125 
- 8.100 
10. 125 
8.100 
- 8.100 



w 

A 

w 

A 

— 

a i p i 

a 2 P 2 

-729 

0 

3.8549 

0 

2025 

0 

-6. 4249 

0 

-2025 

0 

-6.4249 

0 

729 

0 

3.8549 

0 

1296 

1296 

5.1399 

-5. 7810 

0 

-729 

0 

-4. 3357 

0 

2025 

0 

7.2262 

0 

-2025 

0 

7.2262 

0 

729 

0 

-4.3357 

0 

-1296 

0 

-5. 7810 


^ign chosen to give positive 


1.1215 

.7500 


(^f 56 ©*®©x© 



2.1214 

.1515 

.2188 

1.8771 

.0005 

.4729 

.0599 

.0038 

.0139 


m: 



-314. 3 
-134.7 


-.4 

957-7 

-121.3 

2.8 

-18.0 


820.8 

351.8 


a(D = 

a (o) + r A ]-i 134.7 


- 351.8 

a l<D _ 

- 6.4249 0.2707 -6.1542 


+ = 

a 2 (D 

7.2262 -0.4588 6.7674 
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Member 


a 

b 

c 

d 

e 

f 

g 

h 

i 

J 


TABLE IV. - COMPUTATIONS FOR FINAL APPROXIMATION 




Cl } 


- To ) 



CD 

a l p l 

V2/ 

a 2 P 2 

P 

P R 

t p \6. 56 

©*© 


3.7855 

0 

1.1146 

2.0378 

-1485.5 

0 

-6. 3092 

0 

• 7643 

.1715 

347.2 

0 

-6.3092 

0 

.7789 

.1942 

-393.2 

0 

3.7855 

0 

1. 1076 

1.9550 

1425.2 

0 

5.04-74 

• 

-5-5489 

.0619 




0 

-4,-1617 

.3297 

.0007 

0 

-.5 

0 

6.9362 

.8563 

.3615 

0 

732.0 

0 

6.9362 

.6869 

.0851 

0 

-172.4 

0 

-4.1617 

.4110 

.0029 

0 

2.1 

0 

-5-5489 

.5495 

.0197 

0 

-25.5 


a 1 ( r “ 1 ) = -6.3092 
a2^ r “ 1 ) = 6.9362 


I 


MI 


-106.3 


-45.56 


535- T 


229.6 


aW 

= 

-6.4249 


45.56 



7. 2262 

-229.6 


a i (r) 

*2 (r> 

= 

-6.4249 

+ 

0.1157 


-6. 3092 



7. 2262 


-0.2900 


6.9362 



















